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Abstract. We give a new representation of fractional Brownian 
motion with Hurst parameter H < ^ using stochastic partial differ- 
ential equations. This representation allows us to use the Markov 
property and time reversal, tools which are not usually available 
for fractional Brownian motion. We then give simple proofs that 
»^ ' fractional Brownian motion does not hit points in the critical di- 

, mension, and that it does not have double points in the critical 

dimension. These facts were already known, but our proofs are 
quite simple and use some ideas of Levy. 

> 

'nJ- ' 1. Introduction 

m 

^ ■ Our main result is a new representation for fractional Brownian mo- 

^ . tion using stochastic partial differential equations, described in this sec- 

^ ! tion. As an application, in Section 2, Theorems [T] and [2], we state some 

<yQ [ known results about when fractional Brownian motion hits points and 

O ■ has double points. Our representation allows us to give simple proofs 

of these results. 

^ i In recent years there has been an upsurge of interest in fractional 

^ I Brownian motion, see Nualart, Chapter 5 |Nua06] . The most common 

model for noise in physical systems is white noise Bf, the derivative of 
Brownian motion. The central limit theorem gives some justification 
for using a Gaussian process such as Bt. Furthermore, Bg, Bt are inde- 
pendent if s 7^ t. In many situations, however, there are correlations 
between noise at different times. A natural correlated Gaussian model 
to consider is fractional Brownian motion Xt = : t > taking 
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values in R"^, with Hurst parameter H G (0, 1]. The process Xt is 
uniquely specified by the following axioms. 

(1) Xq = with probability 1. 

(2) : t > is a Gaussian process with stationary increments. 
That is, for t,h > the probability distribution of the increment 
Xt+h — Xt is independent of t. 

(3) For c > we have X^ = c^Xt, where = denotes equality in 
distribution. 

(4) Xi has the standard normal distribution in R". 

Note that Brownian motion is a fractional Brownian motion with Hurst 
parameter H = 1/2. 

Next, we describe a seemingly unrelated process, the solution of the 
heat equation with additive Gaussian noise. Then we show that frac- 
tional Brownian motion can be recovered from this solution. There 
are several representations of fractional Brownian motion, see Nualart 
|Nua06j . Chapter 5. One advantage of our representation is that we 
can use the Markov property and time reversal, tools which fail for the 
fractional Brownian motion alone. Using these extra tools, we give a 
simple proof of some hitting properties of fractional Brownian motion, 
and a result of Talagrand about double points. Throughout the paper 
we will write SPDE for "stochastic partial differential equation" . 

Informally speaking, we wish to consider solutions u{t, x) : (t, x) G 
R X R" to the following equation. 



dfU = Au + Fit, x) 
u{—oo,x) = 

where F{t,x) is a generalized Gaussian field with the following covari- 



ance: 



E F{t,x)F{ 



s,y) = 5{t - s)h{x - y) 



where 




and 



H 



1 — an 



2 



so that 



1-2H 



a = 



n 



We will require 



0<H<1 
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SO that 

< a < - 

n 

Our goal is to show that Xt = u(t, 0), but this will not be literally true. 

The above description is not rigorous. To be precise, F is a centered 
Gaussian random linear functional on C^(R""'"^), the set of infinitely 
differentiable functions with compact support on (t, x) G R x R", with 
covariance 

(1.1) Qif,g):=E[Fif)F{g)]= [ [ /{t, x)g{t, x)dxdt. 

ifH=^, and 
(1.2) 

Qif, g) := E [F{f)F{g)] = [ [ [ f{t, x)g{t, y)h{x - y)dydxdt. 

JR" JR" 

if if 7^ Note that in either case, the integral in (11.11) or (11.21) 

is nonnegative definite. Thus, we can extend F{f) to all functions / 
satisfying 

Q(/,/)<oo. 

We call this class of functions H. Note that H implicitly depends on 
a, n. 

Next, for t > and x G R" let 

G(,,):=/(«)-"'^-p(-^) '">0 
[ if t < 

be the heat kernel on R"'. 

We would be tempted to define u(t,x) by 

u{t,x)= / / G{t — s, X — y)F{dyds) 

J-oo JR" 

but the integral will not converge. However, u{t, x) —u{0, 0) looks more 
promising. When if = 1/4 the stationary pinned string was defined in 
pT02j as 

(1.3) U{t,x):=[ [ [G{t-s,x-y)-G{-s,-y)]F{dyds) 

J -oo JR" 

This definition also works for other values of H, provided g eH, where 
9{s, y) = gt,x{s, y) ■= G{t -y,x-y)- G{-s, -y). 

Lemma 1. Let g he as in the previous paragraph. For all t > we 

have 

gis,y)i{s<t) en 
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We will prove Lemma [H in the Appendix. 

From the covariance of F one can easily deduce the following scaling 
property. We leave the proof to the reader. 

Lemma 2. The noise F obeys the following scaling relation. 



F{ct, c^/^x 



V 



'(l+an)/2 



F{t,x) 



Turning to the SPDE, define v{t,x) by av{t,x) = u{ct, c^^'^x). The 
reader can verify the following calculation using fll.3l) . 

adfV = cdtU 

= c(^AU + F{ct,c^/'^x) 

where the equality in distribution holds for the entire random field 
indexed by t, x. Thus, we can cancel out the constants c, a provided 



a = c 2 

and then v satisfies the same equation as u. Thus, 

aU{t,x) = U{ct,c^^^x) 
Setting X = gives us the scaling relation for U{t,0). Thus we find 

Lemma 3. U{t,0) obeys the following scaling relation. For c > we 
have 

U{ct,0) = c^U{t,0) 

where the equality in distribution holds for the entire process indexed 
byt. 

Remark 1. Let Vt^x{s,y) = U(t + s,x + y) — U{t,x). It follows im- 
mediately from M.3\) that the random fields Vt^x{s-,y) and U{s,y) are 
equal in distribution. 



Let 
where 

and 

We claim that 



Xt = K^U{t,0). 

1/2 



(2-a)r( 



if a ^ 1 



Ko. = 2-^/\4nf^ if a = L 
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Proposition 1. Assume that < a < ^, and let Xt = U{t,0). Then 
Xf, as defined above, is fractional Brownian motion with Hurst param- 
eter 




— an 



Proof. We only need to verify the four axioms for fractional Brownian 
motion. It follows from (11. 3p that Xq = 0, so axiom 1 is satisfied. 
Axiom 2 follows from Remark [1] Axiom 3 follows from the scaling 
properties of fractional Brownian motion and Lemma[3l Finally, Axiom 
4 follows from (11. Sp and the integral of the covariance h, which we verify 
in the Appendix. I 



2. Critical dimension for hitting points, and for double 

POINTS 

The rest of the paper is devoted to the following questions. 

(1) For which values of d, H does Xt hit points? 

(2) For which values of d, H does Xt have double points? 

Recall that we say Xt hits points if for each z G R", there is a positive 
probability that Xt = z for some t > 0. We say that Xt has double 
points if there is a positive probability that Xg = Xt for some positive 
times t s. Here are our main results. 

Theorem 1. Assume < H < ^, and that jj is an integer. For the 
critical dimension n = jj, fractional Brownian motion does not hit 
points. 

Theorem 2. Assume < H < ^, and that ^ is an integer. For the 
critical dimension n = jj, fractional Brownian motion does not have 
double points. 

In fact, Talagrand answered the question of double points in |Tal98j . 
Theorem 1.1, and the assertion about hitting points was already known. 
Techniques from Gaussian processes, such as Theorem 22.1 of |GH80j . 
can usually answer such questions except in the critical case, which is 
much more delicate. The critical case is the set of parameters d, H 
which lie on the boundary of the parameter set where the property 
occurs. For example, d = 2, H = | falls in the the critical case for 
fractional Brownian motion to hit points. But for H = ^ we just have 
standard Brownian motion, which does not hit points in R^. This 
illustrates the usual situation, that hitting does not occur, or double 
points do not occur, in the critical case. 
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It is not hard to guess the critical parameter set for fractional Brow- 
nian motion hitting points or having double points. Heuristically, the 

range of a process with scaling Xat = a^Xt should have Hausdorff 
dimension if Xt takes values in a space of dimension at least 

For example, Brownian motion satisfies B^t = o^^^^Bf, and Brownian 
motion has range of Hausdorff dimension 2, at least if the Brownian 
motion takes values in R" with n > 2. The critical parameter of H 
for a process to hit points should be when the dimension of the range 
equals the dimension of the space. Thus, the critical case for fractional 
Brownian motion taking values in should be when the Hurst pa- 
rameter is H = 1/n. For double points, we consider the 2-parameter 
process V{s,t) = Xt — Xg. This process hits zero at double points of 
Xt, except when t = s. The Hausdorff dimension of the range of V 
should he n = j^, and so the critical Hurst parameter for double points 
of Xt should be if = 

First note that the supercritical case can be reduced to the critical 
case. That is, if Xt = (Xt^^), . . . , x^*-""^™-*) = 0, then it is also true that 
the projection (X^^-*), . . . ,Xt"'^) = 0. Furthermore, the subcritical case 
is easier than the critical case, and it can be analyzed using Theorem 
22.1 of Geman and Horowitz |GH80] . Therefore, we concentrate on the 
critical case H = 1/n. 

Below we give a simple argument inspired by [MT02j which settles 
the critical case. The argument goes back to Levy, and an excellent 
exposition is given in Khoshnevisan |Kho03j . It is based on scaling 
properties of the process, the Markov property, and time reversal. Al- 
though fractional Brownian motion is not a Markov process, U{t, x) 
does have the Markov property with respect to time. Furthermore, it 
is time-reversible. 



3. Summary of Levy's argument 

Here is a brief summary of Levy's argument that 2-dimensional 
Brownian motion does not hit points. Let m{dx) denote Lebesgue 
measure on R" and let Bt denote Brownian motion on R". For this 
section, let n = 2. Furthermore, let B[a,b] := {Bt : a < t < b}. It 
suffices to show that 



(3.1^ 



E 



m(E[0,2]) 
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since then we would have 

^ E 



I l{z e B[Q,2\)dz 
P{z e B[Q,2])dz 



and so P{z G i?[0, 2]) = for almost every z. 
Next, for < t < 1, let 

Yt = B,+t-B^ 
Zt — Bi^t — Bi 

Recall that Yt, Zt : < t < 1 are independent standard 2-dimensional 
Brownian motions. This is a standard property of Brownian motion, 
which can be verfied by examining the covariances of Yt, Zt : < t < 1. 
Then F[0, 1], Z[0, 1] are independent random sets. Furthermore, by 
Brownian scaling and translation 



E 



m{B[0,2]) = 2E m(S[0,l]) 



^2E 



m{B[l,2]) 



E 


m{B[0,l]) 


+ E 


m(5[l,2])j 


E 


m{Y[0, 1]) 


+ E 


m{Z[0, 1]) 



On the other hand, set theory gives us 



E 



m(S[0,2]) ^E[m{Y[0,l]UZ[0,l])] 



E 



m(F[0,l]) 



+ E 



m(Z[0,l]) 



-E 



and therefore 



E 



m(F[0,l] nZ[0,l]) 

By Fubini's theorem, 

(3.2) = E[m(F[0,l] nZ[0,l]) 

^ E ly[o,i](^)lz[o,i](^)^^^ 
= / E[lY[Q,i]{z)lz[o,i]{z)]dz 



m(F[0, 1] nZ[0,l]) 



0. 
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By the independence of Y[0, 1] and Z[0, 1] and the Cauchy-Schwarz 
inequahty, we have 



(3.3) 







> 



E [lY[o,i]iz)lz[o,i]iz)] dz 

R2 

E [1yio,i]{z)] E [lz[o,i]{z)] dz 

{E [ly[o,l]{z)]Y dz 

E [ly[o,i](2;)] dz 



R2 



R2 



(^Em{Y[0, 1]) 



Therefore i?[m(y [0, 1])] = and (13. ip follows from the definition of Y. 

4. Proof of Theorems [T] and [2] 
Now we use Levy's argument to prove our main theorems. 

4.1. Hitting points, Theorem [1], The argument exactly follows that 
in Section [3], except that is replaced by R*^, and by axiom (3), 

for c > 0. However, since Xt takes values in R" and H = 1/n, we still 
have 



E 



m(X[0,2]) 



2E 



2E 



m{X[l,2]) 



m(X[0,l]) 

Recall that m(-) denotes Lebesgue measure in R*^. As before, let 

Yt = Xi^t — Xi 
Zt = Xi^t — Xi 

It is no longer true that Y[0, 1], Z[0, 1] are independent. Now we use 
the fact that Xt is equal in distribution to u{t,0), where u{t,x) is the 
stationary pinned string. Changing the probability space if necessary, 
let us write Xt = u(t,0), and let Ht denote the cr-field generated by 
u{t, x) : X E R". Then we have the following lemma. 

Lemma 4. Let us use the above notation. Then Y[0,1], Z[0,1] are 
conditionally independent and identically distributed given Hi . 



Proof of Lemma^ The lemma is proved in [ MT02] . Corollary 1, for 
if = 1/4, and the proof for other values of H uses similar ideas. 
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To show that Y[0, 1], Z[0, 1] are identically distributed, we merely use 
the definition of Y, Z and make a change of variable. Below, equality 
in distribution means that the processes indexed by t are equal in 
distribution. 



Y 



Xi+t — Xi 

f/(l + t,0) -f/(l,0) 



1+t 



[Gil+t 



y)-G{l-s,-y)]F{dyds) 



R" 



[G{t -s,x-y)- G{-s, -y)] F{dyds). 



-oo JR" 

But by the definition of Zt, 

Zt = Xi^t — Xi 

= f/(l -t,0) - f/(l,0) 



[^(1 -s,x-y)-G{l~t-s, -y)] F{dyds) 



-oo JR" 



V 



[ [ [G{t-s,x-y)-G{-s,-y)]F{dyds) 



and therefore Yt, Zt are identically distributed processes. 

Next we discuss the conditional independence of F[0, 1], Z[0, 1] given 
Til. First we claim that the stationary pinned string f/(t, x) enjoys the 
Markov property with respect to t. This is a general fact about stochas- 
tic evolution equations, and we refer the reader to |Wal86] . Chapter 3. 
It follows that Z[0, 1] is conditionally independent of F[0, 1] given Tii. 

This proves Lemma HI I 

From here we duplicate the argument in Section [3l replacing expec- 
tation by conditional expectation given Tii. Briefly, it suffices to show 
that 



(4.1) 



E 



m(F[0, 1]) 



0. 



But, following the same argument as before, we conclude that with 
probability one. 



E 



m(F[0,l] nZ[0,l]) 



0. 



We leave it to the reader to verify that (13.21) and (13. 3p still hold, pro- 
vided expectation is replaced by conditional expectation given TYi. This 
verifles (14.11) , and flnishes the proof of Theorem [H 
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4.2. Double points, Theorem [21 To show that Xt does not have 
double points in the critical case, we use the same argument, but ap- 
plied to the two-parameter process 

V{s,t) := X{t)-X{s). 

We need to show that V{s,t) has no zeros except if s = t. To simplify 
the argument, we will show that V{s,t) has no zeros for {s,t) G TZ, 
where 

7^:= [0,2] X [4,6]. 

The same argument would apply to any other rectangle whose inter- 
section with the diagonal has measure 0. Let us subdivide 71 into 
4 subrectangles TZi : i = 1,...,4 each of which is a translation of 
[0, 1]^. Again we argue as in Section [3l By scaling, we see that for each 

E [m{V{n))] = AE [m(F(7^^))] . 

Next, let TCi be the cx-field generated by {u{l,x) : x G R"}, and 
suppose we have labeled the TZi such that TZi = [0, 1] x [4, 5] and 
7^2 = [1, 2] X [4, 5]. Thus, as before, for each pair i 7^ j G {1, . . . , 4} we 
have 



E 



Hi 



0. 



m{v{ni)nv{nj)) 

Now in |MT02j . Corollary 1, it was shown that for if = 1/4, l^(7^l), 1^(7^2) 
are conditionally i.i.d. given Tii. For other values of H, the argument 
is very similar to the proof of Lemma HI and we leave the details to the 
reader. 

Therefore, as in Section [3], we conclude that with probability one. 



E 



Hi 



0. 



Also as in Section [3l this finishes the proof of Theorem [2l 

Appendix A. Appendix 

Proof of LemmaUl In the case of if = |, where h{x) = 6{x) is not 
difficult to show that 

Gr{x - y^dydr = 2{A7t)-''/H'^^ 

'0 -'R" 

and we leave it to the reader. 

Moving on to the case of H ^ we calculate the following integral. 

Gr{x — z)Gr{x — z')\z — z'\~"'dzdz' 

R" ^R" 
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We make use the Fourier transform, with the following definition. 



/(O = / f{x)e'''-^dx 

for x,^ e R". 

If f{x — z) = Gr{x — z), g{x — z') = Gr{x — z'), h{z — z') = \z — z'\ 
H{z) = g * h{x — z) = J{x — z), then 

H{z) = [ J{x-z)e"< = e''<J{-^) 



therefore. 



Gr{x — z)Gr{x — z')\z — z'\ "dzdz 



f{x — z){ / g{x — z')h{z — z')dz')dz 

R" JR" 

= / f{x-z)H{z)dz=U*H){x) 
Jr" 

It is known that 

(A.2) m = 9i-0 = e-^l^l' 

and by |Str03j . page 51, we get 
(A.3) 

h{^) = = c(a)|er("""^ where c(a) = - ^ ^ - 



r(f; 



by fOO) and (OOj) . (00) becomes 

Gr{x — z)Gr{x — z')\z — z'\~°'dzdz' 



R" JR" 

To evaluate the integral in R*^, we use polar coordinates. Then |,^| = m, 
d{,i) = u"-~^du, where < m < oo, and thus if a = 2{^Y /V{\) is the 
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area of the unit sphere in R", we get 



cr I e u — 

u 



POO 

-| a_ a_ I -1 

cr2 2r 2 / e *^y2 '-dy 
Jo 



(A.4) = r-r^r ^ 

r(f) 



Integrating (1A.4|) over r from to t, we have 



III Gr{x — z)Gr{x — z)\z — z\ "'dzdz'dr 
Jo Jr" Jh"^ 



^ ^ ^ t~ 

(2-«)r(t) 
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